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Abstract 

Given a collection S of subsets of sonic set U, and M C U, the set cover problem is to find 
the smallest subcollection C C S such that M is a subset of the union of the sets in C. While 
the general problem is NP-hard to solve, even approximately, here we consider some geometric 
special cases, where usually U = di^. Extending prior results ,.B_GiI5- , we show that approximation 
algorithms with provable performance exist, under a certain general condition: that for a random 
subset R C S and function /(), there is a decomposition of the complement U \ Uy^rY into 
an expected /(|i?|) regions, each region of a particular simple form. We show that under this 
condition, a cover of size 0(f{\C\)) can be found. Our proof involves the generalization of shallow 
cuttings |Mat92| to more general geometric situations. We obtain constant-factor approximation 
algorithms for covering by unit cubes in 5R'^, for guarding a one-dimensional terrain, and for 
covering by similar-sized fat triangles in 3?^. We also obtain improved approximation guarantees 
for fat triangles, of arbitrary size, and for a class of fat objects. 

1 Introduction 

Given a collection S of subsets of some set U, and M C U, the set cover problem is to hnd the 
smallest subcollection C C S such that M is a subset of union of the sets in C. In the geometric 
setting, almost always U = For example, M could be a finite set of points, and S a given finite 
set of balls. The family S can be specified implicitly; an example is when S is the set of all unit 
balls. Another interesting example is when M is the set of points in a simple polygon in and S 
is the set of visibility regions of the vertices of the polygon. 

The general set cover problem is hard to solve, even approximately, and the simple greedy 
algorithm has performance very close to best possible for a polynomial-time algorithm, assuming 
a certain widely believed complexity theoretic assumption. |Fei98l ILY94j Even in the geometric 
setting, most versions of the problem are believed to be NP-hard, and indeed NP-hardness has 
been shown for several versions. (In some cases, hardness of approximation has been shown as 
well.) The focus of current work is therefore on obtaining approximation algorithms that run in 
polynomial time. Often one obtains a polynomial-time algorithm guaranteeing a logarithmic factor 
approximation by reducing the geometric set cover problem to the combinatorial set cover problem 
|Chv79irOr74llLov75j . 

In many cases, the approximation factor can be made O(logc), where c is the size of the 
optimal solution. Such a result was achieved for the case of polytope approximation in general 
dimension |Cla93j . by applying the iterative reweighting approach jLit87[ IWelSSI ICla95j to an asso- 
ciated set cover problem. (The reduction of polytope approximation to set cover was observed by 
Mitcheh and Suri |MS95, .) 

Bronnimann and Goodrich |B(T95j showed that a very similar algorithm applies in the general 
setting of set systems with finite VC dimension. |BG95] A key observation of theirs was a connection 
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with e-nets. Consider the subset Ue C U, comprising those points of U contained not just in one 
set in 5, but in at least e\S\ of them. An e-net is a cover for such heavily covered points. (That 
is, the set cover problem is to find the smallest possible l/|5|-net for M = U.) Suppose that the 
family S has a 1/r net of size g{r), for every r with 1 < r < IS"]. The algorithm of Bronnimann and 
Goodrich guarantees an approximation factor of 0{g{c)/c), where c is the size |C| of the optimal 
solution. For many cases where g{r) = 0(r log r) |Cla87l lHW87j . their algorithm gives an O(logc) 
approximation. Moreover, if g{r) = 0(r), such as when 5 is a family of disks in or halfspaces 
in |MS W90l IMat92] . they obtain an 0(1) approximation algorithm. 

There have been a few other interesting instances where the O(logc) factor has been improved 
upon. Some recent ones include an 0{\/log n) approximation factor for covering an isothetic polygon 
(with holes) using a minimum number of rectangles contained in the polygon KROSj . and an 0(1) 
approximation algorithm for guarding an x-monotone polygonal chain |BMKM05 . 

Hochbaum and Maass |HM85j consider the problem of covering a set of points in the plane with 
the smallest number of unit disks. For this and some related problems, they present algorithms, 
that for any e > 0, run in polynomial time and return a (1 + e)-approximation. Since any unit 
disk may be chosen in the cover, the problem has a different flavor from that of covering the points 
using the minimum number of disks chosen from a set of specified unit disks. 

Our Results. We generalize results giving small e-nets for halfspaces |Mat92] to a more general 
setting jnSHn] making a connection between the combinatorial complexity of the union of a set of 
objects and size of a net for the set of objects. Suppose that 5 is a set of objects, say triangles in the 
plane for concreteness. Suppose that there is a bound f{j) > j on the combinatorial complexity of 
the boundary of the union of any j objects from S. (More precisely, we need the number of simple 
regions in a canonical decomposition of the exterior of the union of the j objects to be at most 
f{j).) We then show, in Theorem 12.21 that for there is a 1/r net of size 0{f{r)), for every r < \S\. 
We apply a "repair" or "alteration" technique, using a random sample to divide the problem into 
roughly small subproblems, followed by "repair" step in each subproblem. This approach is similar 
to Matousek's; it has been applied also in a similar way to construct "cuttings." |Cha91j 

As noted, this implies a polynomial time algorithm that guarantees an 0(/(c)/c) approximation 
factor for covering a set M of points using objects from S, where c is the size of the optimal 
cover. jESnH Theorem 3.2] (Note that the result is only interesting for /(r) = 0(r log n); otherwise 
the greedy algorithm could be used.) 

We give several applications of this result. If 5 is a set of fat triangles in the plane, then the 
combinatorial complexity of the union of any j elements of S is 0(j loglogj) |MPS"'"94| . and thus 
we obtain 1/r-nets of size 0(r log log r) for fat triangles. This implies, as stated in Theorem 13. ![ 
a polynomial algorithm for the corresponding set cover problem that guarantees an approximation 
factor of O (log log c). If the triangles in S have roughly the same diameter, then the union of 
any j elements from S has a combinatorial complexity of 0{j) [MPS"*" 94 , and we obtain 1/r-nets 



of size 0(r) and an algorithm for the corresponding set cover problem that guarantees an 0(1) 
approximation. There are other applications in this vein. 

Such cover problems are related to wireless network planning, where the sets in S correspond 
to antenna coverage areas. Prior work has sometimes approximated the coverage areas as circular 
disks. |nMWZ04j but often such an idealized model would be far from ideal. Thus the results for 
more general "fat' objects reported here are relevant. 

Another problem that can be viewed as a special case of wireless network planning is that of 
guarding a one-dimensional terrain. Here, the problem is to guard the region above an x-monotone 
polygonal chain using the minimum number of point guards, who are constrained to be on the 
chain. The problem was recently studied by Ben-Moshe et al. |BMKMn5] who presented a fairly 



2 



sophisticated polynomial time algorithm that guarantees an 0(1) approximation. We show that a 
different polynomial-time constant-factor approximation algorithm can be derived quite naturally 
from our paradigm. The approximation result is Theorem 13.51 and applies a generalization of the 
"Order Claim" of |BMKM05] to show, in Lemma 13.41 that an associated sequence is Davenport- 
Schinzel. 

We next consider the case where is a set of axis-parallel unit cubes in 3f?^. Boissonat et 
al. |BSTY98] have shown that the combinatorial complexity of the union of j such cubes is 0(j). 
Such a bound is however not readily available for a canonical decomposition of the exterior of the 
union. We nevertheless exploit the fact that all the cubes have roughly the same size to obtain a 
1/r-net of size 0{r) and, as stated in Theorem 13.81 a polynomial algorithm for the corresponding 
set cover problem that guarantees a factor of 0(1). 

2 General results 

2.1 Small e-nets from small 0-region sets 

In a geometric setting, the set cover and e-net problems often have the helpful structure that for 
any collection H C S, the complement i]\L){H) has a canonical decomposition into locally defined 
pieces. (Here D{H) is short-hand for Uj^g/zy.) That is, there is a set J^{S) of subsets of U, such 
that for any H C S, U \ U{H) can be expressed as a union of sets U \ D{H) = [j{To{H)), where 
J-q{H) C J'{S). Moreover, there is some integer b so that such decompositions J-q{H) can be 
described as follows: for each y G ^{S), there is a configuration By C S of size at most b, such that 
y £ J-q{H) only if By C H and y n D{H) is empty. Say that By defines y in that case. If y n s is 
not empty, for some s G S, say (as usual) that s meets y. So y is in !Fq[H) only if no s € 5 meets 

y- 

It sometimes happens that for some y G ^{S) there is more than one natural configuration By 
that defines y. To reduce problems with such degenerate situations, it is often helpful to consider 
the regions not only as subsets of U, but as configurations {y,By), where By defines y. Also, the 
condition that s £ S meets y will have an analog for configurations, such that s meets or conflicts 
with {y,By) not only if s n y is nonempty, but also if s takes precedence over a member of By, 
for tie-breaking or other reasons specific to an application. The set Tq{H) will be generalized to 
comprise such configurations, and a configuration {y,By) S Tq[H) if and only if By G H and no 
s € H conflicts with {y,By), in this broader way. Even with this generalization, however, we will 
have l]\U{H) C U(jro(/i)), where here U(.7^o(^)) := ^{y,B )<^To{H))y- We may confuse {y,By) with 
y at times, but the situation should be clear in context. 

We will call the configurations in Tq{H) 0-regions. The "0" in Tq{H) and in 0-region indicates 
that the regions do not conflict with the objects in H. More generally, there could be y G ^{3) that 
have By C R, but {y,By) conflicts with j members of H. In that case, say that (y,By) E Tj{R), 
that is, {y,By) is a j-region of R. Note that {y,By) might be a 0-region with respect to R, but a 
j-region with respect to S, that is, conflict with j members of S. 

Call a given combination of U, objects S, regions ^{S), parameter b, defining relation, and 
conflict relation a configuration system. We are assuming that any point not in i? C 5 is in some 
0-region of R. In such a case, say that the configuration system is complete. 

This decomposition of the complement puts the problem into the "object /region" framework [CS89| 
ICMS93j . which is similar to the starter/stopper framework of Mulmulev. |Mul93] Several properties 
of the problem follow from that framework. A basic property within the framework is the following 
version of e-nets, proven in the objects/regions framework. |Cla87] and also in the framework of 
bounded VC dimension |IIW87j . 
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Lemma 2.1 (Likely e-nets) For a given complete configuration system, there is a constant K such 
that, for a random subset R <Z S of size Krlogr, with probability at least 1 — 1/r, every 0-region 
of R is a {<n/r)-region with respect to S, that is, a j-region with respect to S for some j < n/r. 

Since our assumption here is that a point not covered by R is in some 0-region of R, it follows 
that R satisfying the condition of the lemma is an e-net, for e = 1/r and < Krlogr. Call an 
e-net under such conditions a likely e-net. (See Section ^ for the definition of an e-net.) Note that 
by repeatedly sampling an expected 1 + 0(l/r) times, a likely 1/r-net can be found; also note that 
an algorithm for verifying the e-net condition would be needed to apply the lemma. 

Proof: See |Cla87j : also, since the region here have finite VC-dimension, the similar results of 
|HW87j apply. The proof is simply the union bound, applied to every {y,By); the probability is 
small that a particular j-region of S, with j > n/r, is a 0-region of R, and there are 0{n^) j-regions. 
■ 

We will need the existence of such likely e-nets under slightly stronger conditions, which are 
most conveniently stated simply by requiring that they exist for any subset of S. 

Using the existence of likely e-nets, and the objects /regions framework, we have the following 
generalization of "shallow cuttings". 

Theorem 2.2 For a given complete configuration system, let f{r) := E\Tq{R)\, where R <Z S is a 
random subset of size r. Suppose that likely e-nets exist for any subset of S. Then given r > 2b, 
there is a 1/r-net of size 0{f(r)). 

Proof: The construction is as follows. Pick a random subset R' C S of size r. For each 
y G J^o{R'), suppose y meets a set S' C S, of size j'n/r. If / < 1, let let Ry := 0; other, let Ry be 
a likely (l/j')-net for the objects conflicting with y. Such an Ry will have size at most Kj' log j'. 
Then R := R' U \Jy^jr^^(^pr^Ry is a 1/r-net for S, by construction. The expected size of R can be 
bounded using Theorem 3.6 of ,CS89 with c = 2, and the "work" of that theorem is Kj' log j' for a 
( j'n/r )-region, or no more than W{{/!^), where j = j'n/r, and Wi) is the concave "work" function 

2 

W{x) := 4i^^y^log(x^), giving a bound 

^(^((7^^2,5))/(r) = 0(/(r)), 

assuming b is constant, implying also that the term i^2,6 of the theorem is constant. ■ 
We note that the proof suggests a natural randomized algorithm to compute a net. Under 
appropriate assumptions that certainly hold for the applications in this paper, the expected running 
time of this algorithm is polynomial in the input size. 

2.2 Small covers from small e-nets 

Theorem 2.3 For a given complete configuration system, with f{r) as in the last theorem, suppose 
there is a cover C C S of size c for subset Met/. Then a cover o/M of size 0{f{c)) can be found 
in the time proportional to that needed to construct an 0{\/ c)-net, as in the last theorem, times a 
polynomial in \S\. 

(Note that for particular instances a stronger time bound can be obtained.) 

Proof: The previous theorem implies the existence of 1/r-nets of size 0{f{r)). This theorem 
then follows from Theorem 3.2 of |B(T95j . In the algorithm given to prove their theorem, e-nets 
are found many times, for slightly different sets. An alternative approach is to solve the linear 
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programming relaxation, and find a single e-net. [ERS] One version of the latter approach is roughly 
as follows: solve the linear programming relaxation of the problem, which yields an assignment, 
for each object in s G S", of a value Ws with < < 1, such that for each point p £ M, it holds 
that '^p^gWs > 1. Then create a multiset S' , with "copies" of each s, where the number of copies 
is proportional to Ws- Extend the conflict relation with tie-breaking to allow at most one copy to 
contribute to the definition of a region. The resulting configuration system has the property that 
every point in M is contained in \S'\/c regions; that is, a 1/c-net is a cover. ■ 

3 Applications 

3.1 Covering by Fat Triangles or Regions 

Our first applications of the general results follow fairly directly from existing combinatorial bounds 
and the low complexity of trapezoidal decompositions in the plane. 

Theorem 3.1 There is a randomized polynomial time algorithm that, given a set Wlofm points in 
3?^ and a set S of n fat triangles that cover M, computes a subset S' C S o/ 0(c log log c) triangles 
that cover M, where c is the size of the smallest subset of S that covers M. 

Proof: (Sketch) It is long known that the union of n fat triangles has combinatorial complexity 
0(n log log n). (See |MPS+94| . which also gives a definition of fatness.) The same bound applies 
to the canonical trapezoidal decomposition of the complement of their union [Mul93j; we can then 
apply Theorem 12.21 with these trapezoids as the regions. Similar remarks apply for fat triangles 
of approximately the same size, relying on the sharper bounds known for the complexity of their 
union IMPS+Mj ■ 

Theorem 3.2 There is a randomized polynomial time algorithm that, given a set Wl of m points 
in and a set S of n {a,f3)-fat objects of approximately the same size that covers M, computes a 
subset 5" C 5 of size 0(As+2(c)) that covers M, where c is the size of the smallest subset of S that 
covers M. Here, s is the maximum number of intersections between the boundaries of two objects 
in S. 

Here As+2(?^) is a very-nearly linear function of n, related to the complexity of Davenport- 
Schinzel sequences. 

Proof: (Sketch) We use a result of Efrat |Efr99j that the combinatorial complexity of the 
boundary of the union of k such fat objects is 0(Acj+2(fc)), and proceed as in the case of triangles. 
We assume that the trapezoidal decomposition can be efficiently computed (in polynomial time). 
■ 

3.2 Guarding a Monotone Polygonal Chain 

Let P be a x-monotone polygonal chain in 3f?^ with n vertices. Let Q := {^i, . . . ,gm} be a set of 
points, which we will call guards, on P. Say that a guard g lying on polygonal chain P sees a point 
p if the line segment gp does not intersect the region in 3f?^ that is strictly below P. 

Consider the set Mp of points in 3f?^ that are on or above P. For g £ Q, let Ws{g) := {p € 
5?^ I g sees p}, the visibility polygon of g, be the set of all points seen by g. The problem of guarding 
P is that of covering the set Mp by a small subset of S := {Vis(5() | g S G}. For S' C S, the 
complement of the region covered by S' is the area between P and the lower envelope of the visibility 
polygons in S'. Each point on the x-axis has some corresponding point on the lower envelope 
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(perhaps at infinity). It will be helpful, for showing the existence of a a low-complexity, locally- 
defined description of the lower envelope, to consider visibility from the left or right separately. It 
will also be helpful to break ties among the guards determining the lower envelope at a given x 
coordinate. 

3.2.1 Complexity of the Lower Envelope 

Say that g sees p from the left if g sees p and x{g) < x{p), where x{p) is the x-coordinate of point p; 
define visibility from the right analogously. For g € Q, let Lvis{g) := {p £ 5 sees p from the left}, 
the set of points that g sees from the left. Let Sl ■= {Lvis{g) \ g S G}. 

Fix some subset TC Q Q. Given an x-coordinate x, say that g £ 7i owns x from the left (relative 
to 7i) if there is y such that g sees (x, y) from the left, and there is no y', g' £ 7i such that {y' , x{g')) 
is lexicographically less than {y,x{g)). We will say that g owns x from the left at {x,y) (relative to 
TC). If some x-coordinate x is owned by no point in TC, say that x has the owner NULL. 

We can now define the ownership diagram of a set of guards TL ^ G, with respect to P; this 
definition is for ownership from the left, but similar definitions and claims apply for ownership from 
the right. The (left) ownership diagram is the partition of the x axis obtained from the connected 
components of each equivalence class of the relation "x and x' have the same owner." Such 
components are intervals (or single points), and so this diagram is a sequence of intervals, each 
with one owner. Call the corresponding sequence of owners, but excluding NULL, the ownership 
sequence for 7i. A key claim for a bound on the length of this sequence is the following, a slight 
generalization of Lemma 2.1 of BMKM05 . 

Lemma 3.3 Suppose a,b G TC C G and x,x' E have x{a) < x(6) < x < x'. Suppose also a owns 
X (relative toTL) at a point p, and p' = {x',y') is seen by b. Then p' is seen by a also. 

Proof: Since a owns x at p, a sees p, and so P is not above line segment op. Since 6 is on P 
and between a and p, b in particular is not above op. Similarly, P is not above segment bp'. Also 
p is not above bp': if p were above bp', it would be seen by b, and since P is not above bp', b would 
also see some point below p, but with the same x coordinate, contradicting the assumption that a 
owns X at p. So b and P are not above op, and p and P are not above bp'. Therefore a sees p', as 
claimed. ■ 

Lemma 3.4 An ownership sequence for a set TC of r guards is an (r, 2) Davenport-Schinzel se- 
quence, and therefore has length at most 2r — 1. It follows that the number of ownership intervals 
is no more than 2r. 

Proof: An (r, 2) Davenport-Schinzel sequence |SA95j is a sequence of r symbols with no successive 
entries identical, and with no subsequence of the form a ... 6 ... a ... 6. Consider a,b £ Tl, and first 
suppose that x(a) < x(6), as in the previous theorem. It may be that a owns intervals before b 
(with smaller x coordinate than x(6)), and it may be that b owns some intervals to its right, but if 
a owns some x-coordinate at point p, strictly to the right of b, then from the previous lemma, any 
point p' with x{p') > x{p) seen by b is also seen by a. Since x(a) < x(6), such a coordinate would 
be owned by a if either a or 6 owns it, and so could not be owned by b. Therefore, there is no 
ownership sequence of the form a . . .b . . . a . . .b. A similar argument works if x(b) < x(a), and thus 
the first claim of the lemma follows. The length bound for such sequences is long-known |S A95j . 
The final claim follows because there is at most one interval with owner NULL; this is the interval 
to the left of all the guards in TL. ■ 
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3.2.2 Guarding in the objects/regions framework 

We will employ Theorem l2.2l to compute a 1/r-net for the set Sl of size 0(r). (Recall that such a net 
is a subset S' C Sl such that any point belonging to greater than {SlI/v sets from Sl also belongs 
to some set in S' .) In order to apply the theorem, we indicate explicitly how the configurations 
and confiicts are defined. There is a configuration corresponding to every interval in the ownership 
diagram for subsets of Q of size at most 3. Consider an interval / in the ownership diagram of 
{a, b, c} C Q, and suppose b owns each x € I, a owns the interval immediately to the left of /, and 
c owns the interval immediately to the right of I^. The set {a, b, c} defines this configuration. (The 
region of this configuration is the set {{x, y)\x £ I, (x, y) Lvis(6)}.) A guard d £ G \ {a, b, c} can 
conflict with this configuration in two ways: 

1. Relative to the set {a, b, c,d}, d rather than b owns some point x' G /. This of course happens 
if d sees some point with x-coordinate x' that lies below the point p at which b owns x' with 
respect to {a, b, c}. Note that this also happens if d sees p and x{d) < x{b). 

2. Relative to the set {a, b, c,d}, b continues to own all points in / but the interval immediately 
to the left of / is owned by d and not a. Because of the way we break ties in defining 
ownership, this is not a pathological situation at all. A conflict also occurs if d owns the 
interval immediately to the right of I in the ownership diagram of {a, b, c, d}. 

With these definitions, observe that the size of J^o{Ti.), for any subset 7i C Q, is exactly equal 
to the number of intervals in the ownership diagram of 7i, which is 0(17^1) by Lemma 13.41 We can 
therefore use the algorithm of Theorem 12. 21 to compute in randomized polynomial time a 1/r net 
for Sl of size 0{r). We define Sr in a manner symmetric to Sl, and note that the union of a l/2r 
net for Sl and a l/2r net for Sr is a 1/r net for S. 

The above arguments are readily adapted to the case where there can be multiple copies of each 
guard. We can therefore use apply Theorem 12.31 and so a set of guards for P can be found with a 
polynomial time algorithm, of size within a constant factor of optimal. 

Theorem 3.5 Let P be a x-monotone polygonal chain in with n vertices. Let Q := {gi, . . . , gm} C 

be guards, such that Mp is seen by Q . Then a subset C C G that also sees Mp, of size within 
0(1) of optimal, can be found in polynomial time. 

3.3 Covering with Cubes 

We now consider the set cover problem where M is a set of m points in and S" is a set of n 
axis-parallel unit cubes in JR'^ that cover M. We first show that any 1 < r < n, there is a 1/r-net 
for S of size 0(r). That is, there is a subset T C S with |T| = 0(r) such that any point that is 
contained in at least n/r cubes from S is also contained in some cube from T. We also present a 
randomized polynomial time algorithm to compute such a 1/r-net. Prom Lemma l2.H it is possible 
to compute a 1/r-net of size 0(r log r) in randomized polynomial time. 
Let G be the vertices of a grid in of side 1/2. That is, 

G ■■= {(^, ^, ^) I are integers }. 

We "assign" each cube C G S" to some point in G that lies in the interior of C. (Note that 
there is always at least one such point.) Let S[p] C 5" denote the set of cubes assigned to the point 

^If a itself owns the interval immediately to the right of /, then such a configuration would be considered by the 
subset {a, b}. 
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p £ G. For each p £ G such that \S[p]\ > where d > is a suitably large constant, we compute 



a ^.^ i g r i | -net T[p] for ^[p] of size Q^M^Ml) using the procedure described below. Let 



Clearly, 



T:= U r[p]. 

peG;|S[p]|>f: 



m<Eo(^^) = o(dr). 



We argue that T is a 1/r-net for S. Let (7 G be any point that is covered by at least n/r cubes 
from S. Consider the cube E of side length 2 that is centered at q. Each cube in S that contains q 
is contained in E, so it must have been assigned to one of the at most d points in GCiE. It follows 
that there is a point p £ GCi E such that S\p] has at least ^ cubes that contain q. Thus T\p\, and 
hence T, will have a cube that contains q. 



A net for a cluster. We now describe a randomized polynomial time algorithm for computing 
a 1/r-net, for any 1 < r < jS'lp]!, for a "cluster" S\p]. The special property of S\p] is that there 
is a point, namely p, that lies in the interior of all the cubes in S[p]. For any non-empty subset 
S' C S[p], we define a canonical trapezoidation of the boundary of the union of the cubes in S' . 
This is obtained by taking, for each face of each cube in S', a canonical trapezoidation of the 
(isothetic polygon corresponding to the) portion of the face that lies on the boundary of the union 
of S'. Let r(5') denote the canonical set of trapezoids thus obtained. 

Proposition 3.6 For any subset S' C S[p], \T{S')\ = 0{\S'\). 

Proof: Boissonat et al. |BSTY98] show that the combinatorial complexity of the boundary of 
the union of cubes in S' is 0(|5'|). The proposition follows because T{S') is linearly bounded by 
the combinatorial complexity of the boundary of the union of 5'. ■ 

We define the "region" fir corresponding to the trapezoid r G r(S") to be the set of all points 
q £ for which r intersects the segment qp in the relative interior of the segment. It is easy to 
see, using the fact that p lies in the interior of all the cubes in S\p], that the regions {fir \t G r(S")} 
partition the exterior of the union of the cubes in S'. The sets that define and conflict with a region 
fir are defined in the standard way: a cube G £ S\p] will conflict with fj,r if G contains a point in 
fir- We can therefore apply Theorem 12.21 to compute a 1/r net for S[p] of size 0(r). 

Putting everything together, we have the following: 

Lemma 3.7 There is a randomized polynomial time algorithm that, given a set S of n axis-parallel 
unit cubes in SR'^, and a parameter 1 < r < n, computes a subset T C S of 0{r) cubes with the 
property that any point that is contained in at least n/r cubes in S is contained in some cube from 
T. 

It is also straightforward to handle the case where there can be multiple copies of each cube. 
Plugging this lemma into the approach of Theorem 12.31 we obtain the following result for the 
corresponding geometric set covering problem. 

Theorem 3.8 There is a randomized polynomial time algorithm that, given a set M C 3fJ^ of m 
points and a a set S of n axis-parallel unit cubes in that cover M, computes a subset T C S of 
0(c) cubes that cover M, where c is the size of the smallest subset of cubes from S that covers M. 
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We remark that the problem of covering a given set of points by the smahest number of axis- 
parallel unit cubes, where we are allowed to pick any axis-parallel unit cube in our cover, admits a 
polynomial time approximation scheme [HM85j . 

4 Conclusion 

It is worth exploring other versions of the geometric set cover problem where better approximation 
guarantees can be obtained via improved bounds on e-nets. Our work also highlights the need for 
a deeper understanding of the connection between bounds on the union and the size of e nets. 

We close with a natural open problem, which is to obtain polynomial-time approximation al- 
gorithms with a sub-logarithmic guarantee for the geometric set cover problem where M is a set of 
m points in and 5" is a set of n unit balls whose union covers M. 

Acknowledgement. The authors appreciate many helpful discussions with Chandra Chekuri. 
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